MONGE AMPERE EQUATION AND BELLMAN OPTIMIZATION 
OF CARLESON EMBEDDING THEOREMS 



VASTLY VASYUNIN AND ALEXANDER VOLBERG 

Abstract. Monge-Ampere equation plays an important part in Analysis. For 
example, it is instrumental in mass transport problems. On the other hand, the 
Bellman function technique appeared recently as a way to consider certain Har- 
monic Analysis problems as the problems of Stochastic Optimal Control. This 
brings us to Bellman PDE, which in stochastic setting is often a Monge-Ampere 
equation or its close relative. We explore the way of solving Monge-Ampere 
equation by a sort of method of characteristics to find the Bellman function 
of certain classical Harmonic Analysis problems, and, therefore, of finding full 
structure of sharp constants and extremal sequences for those problems. 



1. Introduction 

The goal of this article is to show how Monge-Ampere equation allows us to 
solve a certain class of Harmonic Analysis problems. We choose two problems to 
illustrate the method: John-Nirenberg inequality (JNI) and Carleson embeddding 
theorem (GET). 

In these problems we determine Bellman function of the problem (see later). 
Bellman function carries all the information about the problem: sharp constant, 
construction of extremal functions/extremal sequences; and sometimes it also car- 
ries auxiliary information. 

In both JNI and GET Bellman functions were found before: in [SVj for JNI 
and for GET in [M]. However, the way they were found is rather specific for each 
problem and not "from basic principles" so-to-speak. 

We propose here a universal method which fits many problems (we count at least 
seven of them). Among those we chose JNI and GET bacause of their classical 
nature and because the reader can easily compare our approach and our Bellman 
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functions with those in the literature |SVj . [M] . Our approach will differ completely 
from the previous ones, the functions will coincide. 

We borrowed the Monge-Ampere solution by "the method of characteristics" 
from the literature, and the most important inspiration came from Slavin and 
Stokolos [SiSt] . 

Let us notice that the aim of the paper is not only to present a new and unified 
approach to finding Harmonic Analysis Bellman functions. In fact, for GET we 
plan to prove some new results, not just to reprove the old ones due to [M]. 

1.1. Bellman functions. The method of Bellman functions in Harmonic Analysis 
appeared probably in the 80's at the series of papers of Burkholder devoted to sharp 
L^-estimates of martingale transform. It appeared there under certain disguise 
(and Bellman function and Bellman PDE were not even mentioned). The method 
was rediscovered in 1995 in the first preprint version of [NTVlj . It was very much 
developed in [NT] . The method turned out to be useful. Sometimes it solves certain 
important Harmonic Analysis problems that do not have an alternative (classical) 
approach so far, see e.g. [PVoj . There is now certain amount of literature on 
Bellman function technique. We note here an excellent paper [M], which serves as 
our inspiration here, and two small surveys: |NTV2| and |Vo| . 

We fix e G (0, 1) and an interval /. Subintervals of / will be called J . Here is 
what is called in [SV] the Bellman function for the John-Nirenberg inequality: 



BjNi{xi,X2]e) ■.= s\xp{{e'^)j: {4>) ^ = xi, {4?) ^ = X2, 

{<t^)] < {<t^') J < {<t^)] + VJC/}. (1.1) 

Function Bj^r/ is defined in 

^JNI {x = {xi,X2): x\<X2<x\ + e^}. (1.2) 

Now given any / we consider its dyadic lattice of subintervals D = D[I) (it 
consists of left and right halves (sons) of / and then all left and right grandsons 
and et cetera). Let ^ denote a non-negative measure on / without atoms. We want 
to write the Bellman function for the Carleson Embedding Theorem (dyadic). 
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BcET{xi,X2,X3;m,M) := sup-^^ / \(t){t)\'^ dn{t) + ^ \{4>)jfaj: 

)l = X2, |7| (/^W + X] "-^ 



JG-D(/) 

m|J|<^(J) + ^ a£<M|J| VJg £»(/)}. (1.3) 

Functions BcET is defined in 

ET ■= {x = ixi,X2,xs): xj < X2, m < X3 < M} . (1.4) 

We introduce also B(xi , X2, X3). This is what is often cahed the Bellman function 
for the Carleson Embedding Theorem, see [M], |NTj . 



B(xi,X2,X3) := sup — ^ |(0)jPaj: 

{(p)l=Xl, {(t)^)j=X2, "^)=^3, 

5^ a£<|J| VJGi^(/)}. (1.5) 

1&D{J) 

Remark. The reader can easily see that the definitions of all Bellman functions 
BcET, B does not depend on the interval I. 

The formula for the function 'Bjfyfj{xi,X2) was found in [Vj and [SVj . the formula 
for B(xi, X2, X3) was found in [M], the general formula for Bcet{xi,X2, X3;m, M) 
below is new. Let us also notice that formula for B(xi,X2,X3) was found in [M] 
from intricate combinatorial consideration. We tried to explain in the present paper 
a natural way to find this formula from "the basic principles". 

Namely, all formulae will be obtained here by a unified method of solving the 
Monge-Ampere equation. 

Theorem 1. 

Bjm{xi,X2;e) = (l - - (x2 - xf)") e-i+^/^M^^^. 
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Theorem 2. Consider 

X3X2 

[l-2a]'[l-4a(l-X3)] 
where a = a{x) is the unique solution of the cubic equation 



B{xuX2,xs) = „ ^ ' (1-6) 



^2 



Xi 1 — 2a(l — X3) 
1 - 2a 



1 -4a 



1 - 4a(l - X3) 



[1.7) 



X2 

on the interval [0, |] . Then 

B{x)=B{x) xl<X2, 0<X3<1. 

Theorem 3. Consider 

(x^ — m)x2 

B[xi,X2,X3;m,M) = -j- + mx2 , (1.8) 

[1 - 2a{M - m)] [l - 4a(M - X3)] 

where a = a{x) is the unique solution of the cubic equation 



X 



1 

X2 



1 - 2a(M - X3) 
1 - 2a(M - m) 



1 - 4a(M - m) 
1 - 4a(M - X3) 



(1.9) 



on i/ie interval [0, 4^71:;^] . T/ien 

B{x;m, M) xf < X2, m < X3 < M 



BcET{x;m,M) 



mx2 Xi < X2, X3 = m. 



Remark. The cubic equation (jl.9p has actually sometimes two solutions and 
sometimes one solution (for the left hand side in [0,1)). The one (called a(x) 
above) is always in [0, 4(jv/-m) ] • Another one (if exists) is negative. We will call 
it a~ = a~(x). We will see now that it is also responsible for the meaningful and 
interesting extremal problem. 

1.2. Lower Bellman function. If we denote the Bellman function (|ll.ip by Bj-q^^ 
and introduce another Bellman function 

Bjnin{xi,X2,X3;m,M) := inf I / \4>{t)\'^ diJ.{t) + ^ \ {(/)) j\'^aj : 

■'^ JeD{i) 

{(p)j=Xl, {4>'^)j=X2, jJi (/^W + X] «-^)=^3, 

m\J\ < fi{J) + ai<M\J\ VJg £»(/)}, (1.10) 
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Theorem 4. Consider 
B{xi,X2,X3;m,M) 



{X3 - m)x2 



+ mx2 



[1 - 2a(M - m)] [l - 4a(M - X3)] 
where a = a~ (x) is the unique negative solution of the cubic equation 



(1.11) 



X2 



1 - 2a{M - X3) 
1 - 2a{M - m) 



1 - 4a(M - m) 
1 - 4a(M - X3) 



Then 



Bnim(x; m,M) = < 



M) x\ < X2, M 



{M - m)^ <xs<M 

X2 



mx2 X? < X2, m < X3 < M — (M — m)— . 

X2 



Discussion. 

• 1. We are going to provide a detailed proof of Theorem [3] only. We leave for 
the reader to fill out the details of the proof of Theorem SI where concavity 
should be repalced by convexity. 

• 2. An interesting observation follows from the comparison of ()1.6p and (jl.Sh . 
In fact, we can see that the following rescaling relationship holds 

'BcETixi,X2,X3;m,M) = {M - m) B{xi, X2, — ) +mx2- (1.12) 

M — m 

In principle, (jl.l2p should have followed just from the definitions of B, 
BcET above. Let us provide a small explanation. 

First, introduce the notations, by denoting a/ to be the center of dyadic 
interval /, then 



Q := ai + i— , da := 



ai d5ci , 

IdDil) 

where 6^ stands for the delta measure at z as usual. Actually, (jl.l2p means 
that if we are forced to have the uniform estimate from below on how much 
measure we have in any closed Carleson square: 

(^ + a)(Qj) >m|J|,VJGD(/), 

then to maximize the quantity in the definition of Bc£;j'(x; m, M) over 
all admissible //'s and a's (that is to obtain Bc£;r(x; m, M)) one needs to 
keep "boundary measure" to be exactly equal to mdx. In principle, it is 
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easy to believe that always /x should be > mdx. But (I1.12p claims more, 
namely, the equality. It means, that if we want to maximize the outcome 
in the definition of ^cet{x] m, M) we need to keep minimal possible mass 
(namely, mdx) on the boundary, as a reserve so-to-speak, and we need 
to distribute the rest of mass in the form of measure a somewhere inside 
the half-plane. Also (jl.l2p means that this distribution of a mimics and 
actually equals the rescaled distribution of the extremal measure a in the 
case when we maximize B and not '^cet{x\ m, M). 
• 3. Of course (jl.l2p shows that function B is "more equal" among ^CET- 
But interestingly enough, (|1.12p seems not to follow directly from the def- 
initions of functions B and ^cET{x]ra,M) in (jl.Sp and (jll.ip but rather 
from the proofs of Theorems [2] and [3j In fact, it is a priori not clear why 
the optimal distribution of measure for Bc£;r(rE; m, M) should have these 
two features: a) it leaves minimal possible "reserve" on the boundary of 
the circle; b) it distributes the rest by repeating the best distribution for 
B (with correctly rescaled parameters) as (|1.12p shows. 

We shall not prove Theorem[Tl its proof can be found in [V] or |SVj . in the latter 
paper the comparison of dyadic and non-dyadic cases for JNI can be found as well. 
We use this theorem as a "simple" example to illustrate how to find the Bellman 
function by using the general method of solving the Monge- Ampere equation. Then 
we use the same approach in much more involved situations of Theorem O This 
Bellman function is also not quite new, for the case m = 0, M = 1 the function 
was found by Melas in [Mj, where the reader can find the complete proof of the 
theorem for this the most important case. 

However, in [M] there is no explanation how the author was able to guess his very 
non-trivial function. And the main accent of our paper is just the explanation how it 
is possible to guess right, or more correctly, how to solve the corresponding Bellman 
equation, which is a Monge- Ampere equation in the situation under consideration. 
Nevertheless, we will not stop after finding the Bellman function from Theorem [3] 
and provide the complete proof for the sake of completeness and the convenience 
of the reader and because in such more general setting it is new. 

2. How WE PROCEED? 

The consideration of the theorems above can be split to four parts. 
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• I. In the first part one observes that just by definition Bellman functions 
B satisfy a certain concavity condition in their domain of definition and 
boundary conditions on (part of) the boundary of this domain. 

• II. In the second part one considers all function satisfying these concavity 
and boundary conditions. We denote this class by V. And one makes the 
following supposition: as b belongs to V and is the "best" such function, 
it has to satisfy not only the concavity condition but also this concavity 
should become degenerate, i.e., the inequality has to turns into equality 
along some vector field on our domain fi. This brings to the picture the 
Monge-Ampere equation. One solves it using the boundary conditions 
mentioned above. The result is a function B £ V. Function B carries an 
interesting geometric information to be used later. 

These two steps are in fact not necessary for the proof of the results, they are 
needed only to finding a function B, a candidate for a role of the Bellman function 
B. For example, in the excellent paper [M] very complicated Bellman functions 
appear as deus ex machina. As we shall see the analysis of Monge-Ampere equation 
not only supplies us with a candidate, but it helps in the next two steps as well, 
namely, in the prove that the found candidate really is the desired Bellman function. 

• III. The third part consists of proving that > B. In convex domains of 
definition this is usually not difficult. Otherwise it may require a non-trivial 
proof, see |SVj for non-convex 

• IV. The fourth part consists of proving i? < B. This is achieved by pre- 
senting the extremal functions or extremal sequences of functions. In its 
turn such functions are found from the geometric structure of B (mentioned 
above in II). 

3. MONGE-AmPERE EQUATION. 

This section is here for illustrative purposes. Through this section the function 
B is the Bellman function Bjat/ for the John-Nirenberg inequality defined in (II. ID 
and ^} = QjNj is its domain defined in (jl.2p . 

Let us very briefly review Part I for Theorem [TJ 

It is very easy to see that in the dyadic case (i.e. in the case when only dyadic 
intervals J G D{I) are considered in (|l.ip ) the function B has to satisfy concavity 
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condition 

B(x)-i(B(x+) + B(x-)) >0, (3.1) 

for all triples x, x'^ , x~ such that x = ^""""^^ and all three points are from J7. 

In the non-dyadic case (i.e. in the case when all subintervals J C I are considered 
in (jl.ip ) it is not clear a priori why the Bellman function B must satisfy the concav- 
ity condition (|3.ip . And it turns out that for the Bellman function condition (j3.ip 
is not fulfilled. However, it turns out that the Bellman function is locally concave, 
i.e., concave in any convex subdomain of 0,. We shall not discuss here this subtle 
moment, especially because we shall use just the local concavity condition, which 
is of course weaker than the global concavity, and coincides with the latter only for 
convex domains. The discussion of the difference between dyadic and non-dyadic 
Bellman functions for the John-Nirenberg inequality the reader can find in |SV| . 
Here we only add that for the functions smooth enough, the concavity condition 
can be rewritten in the differential form: 



d^B{x) ■= "^"^ ^^"^ < , Vx G , (3.2) 
where by '^x x- we denote the partial derivatives — — - — . 

^ OXiOXj 

We have one obvious (from the definition) boundary condition on the lower 
parabola of the boundary d^l: 

B(xi,x?) = e"S (3.3) 

because only the constant test functions (p correspond the points x with X2 = x^. 
We also have a simple (from definition) homogeneity condition 

B(xi +t,X2 + 2xit + t^) = e*B(xi, xa) . (3.4) 
It follows from the definition of B if to take (f) + t as a test function rather than (p. 
Let us move to the Part II. Put 

V = {v£ C'^{n) : V satisfies i^} . 

Usually any function from V (and moreover, any function satisfying only the 
concavity condition (j3.ip ) supply us with some estimate. But if we looked for 
a sharp estimate we need to choose the minimal possible function v from this 
class, which must be the Bellman function B. This function, "a candidate in the 
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true Bellman function" will be denoted by the usual letter B. For every point 
X £ Q there exists an extremal function (p (or "almost extremal" function (pn, i.e., 
a sequence of functions) realizing the supremum in the definition (jl.ip . The usual 
procedure of using the Bellman function consists in the consecutive application 
of (j3.ip . when splitting the interval I, where a test function is defined. For the 
extremal test function there has to be no lost in such procedure, therefore, for the 
Bellman function the equality has to occur at least for one splitting the point x 
into a pair {x^, For a concave function this means that it is linear at some 

direction. If we have almost extremal functions, i.e., an extremal sequence, then we 
must have "almost equality" in (jS.ip . at least up to the terms of second order. In 
any case this means that the Hessian matrix ()3.2p has to be degenerated. Thus, we 
are looking for the "best" B, on the top of this condition of negativity of Hessian 
we will impose the following degeneration condition: 



vxg173Ggm2\{o|^ ((d2i?)e,e) = ( (^"'^^^^ ^"^"^(^^,(^^1=0. 



(3.5) 

Since the matrix ^B is negatively defined, we conclude from (|3.5p the following 
degeneration condition on the Hessian: 

\Bx2X\ Bx2X2 J 

We underline once more that in the first two steps we can allow ourself not too 
rigorous arguments and various assumptions, because this is not the proof, it is 
heuristic way of finding a candidate in the Bellman function. The rigorous proof 
that the found candidate is indeed the required Bellman function starts from the 
Part HI. 

Claim. There is a simple algorithm to find the family of functions enumerated by 
a parameter 6 {e < 6 < 1) 

B = Bs{x,,X2) := ^ - ^S^-(x2-xj)^ ^x,+^5^~{-2-xl) 

that solves Monge-Ampere equation ()3.6p with boundary condition p.3p and ho- 
mogeneity condition (j3.4p . 
To prove this claim we need 
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3.1. Monge Ampere equation and method of "characteristics". Let v = 

f (xi , .. ., Xfi ) is a smooth function satisfying the following Monge-Ampere equation 
in some domain 



det dFv = det 



0, Vx = (xi, . . . ,Xn) G , (3.7) 

and suppose that this matrix has rank n — 1, i.e., all smaller minors od (Pv are 
non-zero. Then there are functions ti{xi, . . . , x„), i = 0, 1, . . . , n, such that 

v{x) =to + tlXi + t2X2 H h tnXn (3.8) 

and the following n — 1 linear equations hold: 

dto + xidti + X2dt2 H 1- Xn-idtn-^i + = . (3.9) 

Let us explain why this is n — 1 equations and why they are linear. One needs to 



read ()3.9p as follows: we think that, say, ti, . . . , are n — 1 independent variables 
and t„, to are functions of them. Then p.9p can be rewritten as follows 



— +Xi - dti H h ^- hXn-l +Xn^^ ]dtn-l = 0. 

oti oti J \Otn-l Otn-l' 



whence 



dtn dto . , T 



So we get n — 1 equations. 

Remark. In general we can choose any n — 1 variables as independent, of course. 
Since the order of variables is arbitrary, sometimes the first n — 1 is not the most 
convenient choice. 

Now why these are linear equations? We think that ti,. . . ,tn-i is fixed. Then 
the n — 1 equations give us linear relationships in xi, . . . , Xn, so n — 1 hyperplanes. 

Therefore, ()3.9p gives the intersection of n — 1 hyperplanes, so gives us a line. 
We can call it . These lines foliate domain S7 and (j3.8p shows that u is a 

linear function on each such line. 

Let us prove all these propositions. Matrix d'^v annihilates one vector 0(x) at 
every x = (xi, ...,x„) G Q. So we get a vector field 0. Consider its integral curve 
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xi = xiixn)-, ■ ■ ■ iXn-i = Vector Q{x) is a tangent vector to that curve, 

i.e., 

/ x', \ 



e = e. 



^n-l 
1 



(3.10) 



Consider a new function = f (xi(x„ ),..., a;„). Due to (j3.10p its 

second derivative is 



/ x\ \ 



g" = {d'v 



^n.-l 



/ X' \ 



^Jl-1 



1 



+ VxiXi + - ■ ■ + Vx„_iX^_i — VxiXi + - ■ ■ + Vx„_j^X„_i 



(3.11) 



Now, let us also show that Vx^, ■ ■ ■ ,Vx„_-^ are constants on this integral curve. 
Suppose we are standing on the integral curve. The surface Vx^ = ti = const has 
normal {vx^xn • • • , Vxix„)i that is the first row of matrix d?v, which is orthogonal to 
the directional vector G of the integral curve. Hence Q is in the tangent hyperplane 
to the surface Vxi = ti. The same is true for the surfaces Vx^ = U = const, 
i = 2, . . . , n — 1. Intersection of these surfaces gives us our integral curves, because 
is in the intersection of all tangent planes to these surfaces, therefore the curves 
Cii,...,i„_i enumerated by constants are just the integral curves of the 

tangent bundle B. Thus, (|3.11|) can be rewritten as 



dX rir. 



(g - (tlXi H \-tnXn-l)) = 0. 



(3.12) 



We obtain that the second derivative of a function (of x„) in (j3.12p is zero. So 
function is linear in Xn, that is tnXn + tQ, where the constants t„, to depend only 



on the curve Cf, f , , that is 



tn — tn{tl, ■ ■ ■ , tn-l) j ^0 — *o(^l) • • • ) ^n-l) • 

We obtain on Ct, t i 

t-l ) 1 ■'n — 1 



v{xi{Xn), ■■■ , Xn-liXn),Xn) = to + 2^1*1 H 1" Xn-l 
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Since we assumed our vector field to be smooth and its integral curves foliate the 
whole domain, varying parameters ti, . . . ,tn we get (I3.8p . To check (I3.9p take a 
fuh differential in i^M)- Then 

[Uaj^dxiH \-Vx„dxn] = dv = [dto + tidxi^ \-tndXn]+Xidti-\ \-Xndtn- (3.13) 

We are on Ctj^^,,,^t„ and so f^;. = tj, i = 1, . . . , n — 1 as we established already. But 
it is also easy to see that 



on Ct^^,,,^tn~i- In fact, all f^. and all ti are symmetric. We could have chosen 
to represent the integral curve of not as xi = Xi{xn),i = 1, ... ,n — 1 but as 
~ ^i(^i)' J = 2, . . . , n. Now we see that two expression in brackets in ()3.13p are 
equal. Then we obtain xidti + • • • + Xndtn + dt^ = 0, which is desired n — 1 linear 
relationships (|3.9p . 

3.2. The proof of the claim. We can use the previous Section [3.11 But we prefer 
to repeat it for the simple case of only two variables! 

Proof. Equation (13. 6p of course is solved by linear function in xi, X2, but linear 
function cannot satisfy the boundary condition. So we assume that matrix d^B 
is negatively defined, non-zero, and there exists exactly one vector e G R2 \ {g} 
that annihilates it: d^BQ = 0. This vector may depend of course on (xi,X2) and 
we get a vector field Q{x) in $7. Let us see now that Bx^ji = 1,2, are constants 
on the integral curves of this vector field. Locally our integral curves are level sets 
of a certain function s = s{xi,X2). Of course, the function s{x) is not uniquely 
defined, it is defined up to change of variables s <p o s, where (/? is a function of 
one variable. But in any parametrization the vector (— Sx2) ^xi) is a tangent vector 
to the integral curve, i.e., it is collinear with our vector field and thus annihilates 



Notice that both functions Bx^i^xi, X2) instead of s{xi,X2) satisfy the last equation 
as well. In fact, the vector {—{Bxj^)x2j {^xi)xi) annihilates the first row of d^B 
obviously. Then it annihilates the second row just by (j3.6p . The same is true for 
the vector (— (^^.J^j, (i?2,'i)xi)- Therefore, the integral curves of our vector field 
are the level sets of the functions both Bx^{xi, X2) and Bx2{xi,X2)- Thus, we can 




tn-l) 



d^B: 




(3.14) 
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conclude that (at least locally) i?^. = ti o s for some functions ti of one variable. 
Moreover, this means that we can take any of B^^ to play the role of s: to define 
the integral curves of our vector field Q as its level sets. 

Let us check now that the function to := B — tixi — t2X2 is also constant along 
any integral curve. For this aim we have to check that we can put to instead of s 
into ()3.14p . but it is evident, because at each point {xi,X2) we have 

dto dB dti dt2 

-T{ — = - t: \- t: — 3^1 + 7^ — X2 + t2 = -Bxixaa^i + B^2^^X2 

0x2 0x2 0x2 0x2 

dto__dB__dh _dt^ __R _R 

So, we obtained 

B{x) =to + hxi+t2X2. (3.15) 

The integral curve of the kernel field of the Hessian of this function are the straight 
lines given by the equation 

dtQ + xidti + X2dt2 = . (3.16) 

Indeed, on the one hand 

dB = B^-^dxi + Bx^dx2 , 

on the other hand the differential of (j3.15p is 

dB = dtQ + xidti + tidxi + X2dt2 + t2dx2 ■ 

Comparison of the these two equalities yields (j3.16p . 

If the vertical integral curves of our vector field are nowhere dense (so they could 
be given by X2 = X2(xi)), then they could be considered as the level sets ^2 = t, 
and two other parameters are functions of t. Equations (j3.15p and (j3.16p could be 
rewritten as 

B{x) =toit) + tiit)xi+tx2, to{ty + ti{tyxi+X2 = 0. (3.17) 

If the horizontal integral curves of our vector field are nowhere dense, they could 
be written by xi = Xi{x2), or as the level sets ti = t, and two other parameters 
are functions of t. Equations (j3.15p and (|3.16p then takes the form 



B{x) =toit)+txi + t2it)x2, to{ty + xi+t2ityx2 = 0. (3.18) 
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Till now we have considered general equation when the determinant of the Hes- 
sian for a function of two variables is zero. Now we return to our specific problem 
and use (j3.4p to determine the structure of the family of integral lines. 

Differentiating (|3.4p in t and setting t = we get 

or 

ti + 2X1*2 = + tl^l + t2X2 . 

This is the equation of a straight line 

{to - ti) + Xi{ti - 2t2) + X2t2 = G (3.19) 

representing a level set of the functions tj, therefore it has to be the same straight 
line that is given by equation ()3.16p . Hence the following two equations must be 
true 

dto dti dt2 

to — ti ti — 2t2 t2 

They could be rewritten as 
Solution of ((3:211) is 

ti = -2t2log|i2| +2cit2 (3.22) 
(we denoted the integration constant by 2ci for the future convenience). After 
plugging this solution into (|3.20p . we get 

to = t2 log2 |t2| - 2cit2 log \t2\ + C2t2 . (3.23) 

Now we rewrite equation ()3.19p using the obtained expressions for ti: 

log^ \t2\ + 2(1 - Ci) log \t2\ + (C2 - 2ci)) - 2X1 (log |t2| - Cl + l) + X2 = . 

It is convenient introduce a new variable 

a := log|t2| - ci + 1. 
Then the formulas (j3.22p and (j3.23p can be rewritten as follows: 

ti =2(1 - a)t2 , (3.24) 
to={{l-af-cl + C2)t2, (3.25) 
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and the extremal lines are 

X2 - 2axi + + C2 - 1 - c? = . (3.26) 

Since 

c\ + 1 — C2 = X2 — 2axi + = (x2 — x\) + (xi — of' > , 
we can introduce a new non-negative constant 6 such that 

and rewrite equation (j3.26p in the form 

X2 - 2axi + - 5^ = . (3.27) 

We observe that these are lines tangent to parabola X2 = x\ + 5"^ at the point 
xi = a, X2 = a? +6"^ and intersecting the lower boundary X2 = x^ at the points with 
xi = aib(5. We want these lines to foliate domain Q. = {(xi, X2) '■ x\ < X2 < x^+e^}. 
Hence, 

5>e. 

Now we can almost write down B[xi,X2)- The variable a can be found in terms 
of Xi by solving equation (|3.27p : 



a = xi + \l 5"^ - {x2 - x\) . 

We have chosen the solution with the plus sign in front of the square root. The 
explanation, why does it correspond to our problem, and to what extremal problem 
corresponds the opposite sign, can be found in [Vl ISV| . Using now (j3.24p . (j3.25p . 
and (I3.27P we obtain the expression for B: 

B{xi,X2) =to + tlXi + t2X2 

= ((1 - af - c? + C2 + 2(1 - a)xi + X2)t2 

= 2(1 - a + xi)t2 = ±2(1 - a + 2;i)e"+'=^"^ 

= ±2(1 - ^<52 - (3:2 - xf))e^i+V'5'-(^2-^-?)+ci-i ^ 

Since B has to be positive and to satisfy boundary condition (j3.3p we choose the 
sign and the value of the constant ci. Finally, we get a family of solutions 



B{xr,X2,6) = W^!ESlZMIe-+V^M^;=^-^ 

1 — 

depending on a parameter 5. Acceptable values of the parameter 5 are < 5 < 1. 
Why it is so, as well as the choice of this parameter for a given e (by the way. 
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different for the dyadic and non-dyadic cases), an explanation of all these facts can 
be found in |SV] . □ 

The above family of functions B was found in [V] and [SV] from different rea- 
soning. After finding this candidate the steps I, II are completely finished for 
John-Nirenberg inequality. The reader is referred to [SVJ , where steps III and IV 
are done. The fact that Q is not convex makes step III rather delicate. 

Now we are going to present for the reader all four steps for another more difficult 
problem of finding the Bellman functions of Carleson embedding theorems (GET). 

4. Step I for Carleson embedding theorems 

From the definition of function B above we can immediately see that it satisfies 
the inequality (just concavity) 



B(a;) - -(B(x-) +B(a;+)) > 0, Vx=^ G , x = -{x+ + x') . (4.1) 

Note that now i7 = flcET (see (II. 4p ) is convex and for all x^ ft we have x £ fl. 
For the sake of brevity, we shall usually omit the parameters m and M. 
The boundary conditions also follow from the definition: 

B{xi,xl,X3) = xjxs. (4.2) 



— {xi,X2,M;m,M)=xl. (4.3) 
0x3 

It is equally easy to see that homogeneity condition holds 



B{txi,t^X2,X3) =t'^B{xi,X2,X3), tGM. (4.4) 

The detailed explanation of these properties can be found in [NT], where this 
Bellman function was defined and a majorant (so called supersolution) was found. 

Let us consider all smooth functions in 0, satisfying (14. 1|) . ()4.2p . ()4.4|) . Call this 
family V. 

As before, we are looking for the "best" B £ V, so on the top of the con- 
dition (j4.ip of negativity of Hessian we will impose the following degeneration 
condition: 
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Vx = (xi,X2,X3) G 3eeM^\{0} 





/ 




BxiX2 


B 


((^^5)6,9) = 






Bx2X2 


B 




I 




Bx3X2 


B 



XlX'i 



X2X3 



X3X3 




0. (4.5) 



As matrix (PB is negatively defined we conclude from (|3.5p the following degen- 
eration condition on the Hessian: 







^ BxiXl 


Bx\X2 


Bxix:i 








det{d'^B) 


= det 


Bx2Xi 


Bx2X2 


Bx2X3 


1- 




(4.6) 






\BX3XI 


Bx3X2 


BX3X3J 









Claim. There is a simple algorithm to find the function B{x;m, M) that solves 
Monge- Ampere equation (j4.6p in the domain (jl.4p with boundary conditions (|4.2p - 
(j4.3p and homogeneity condition ()4.4p . 



Proof. Again we consider the vector field such that d'^B{x)Q{x) = 0, x = 
(xi,X2,X3) G Vt. And we consider its integral curves. Section l3.ll shows that 
these integral curves are straight lines (segments of straight lines). Section 13.11 
and its the method of characteristics allows us to write the following "parametric" 
equation for these lines, which we first write in the invariant form (compare with 

dSlSD) 

xidti + X2dt2 + xsdts + dtQ = . (4-7) 

This can be rewritten as follows if we choose (ti,t2) as the set of independent 
parameters defining lines -Z^t^.tj foliating our domain: 



dt3 ^ j_ dto _ 

(4.8) 



n -I- HQ. . tq -I- ^ — n 



And similarly to the previous section the solution of the Monge-Ampere equation 
is given by 

B{xi,X2,X3) = ti ■ Xi + t2 ■ X2 + ts ■ xs + to . (4.9) 
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Here each line is given by fixing two free parameters (^1,^2)) wliich are ti = 

. (This can be obtained exactly as in the previous section.) Parameters 



t2 



dB 

8x2 



t3,to are not free, they are unknown functions of ti, t2, for which we will find (a 
priori non- linear) PDE. They will be easy in our case (and linear). And we will 
solve them easily. This will bring us the formula for B in the same way we get 
when proving the claim of Section [3l 

Let us use (|4.4p now. Differentiating in t and setting t = 1 we get 



Whence, 



xiti + 22:2*2 = 25. 



+ 2*3X3 + 2*0 = 0. 



(4.10) 



(4.11) 



Homogeneity gave us (|4.1ip and we conclude that this plane (for each fixed ti, 
t2) contains the line (j4.8p . But the line (j4.8p is passing through the point 



and has the direction 



Hence we get from ()4.1ip : 



. dtp dtp 



, dt3 dts 

^ dti' dt2' 



(4.12) 



(4.13) 



dU 

oti 
dtp 

2to-h-^ = 0. 

Oti 

From these "PDE" we easily write down 

t3 = A{t2)tl 

tp = D{t2)tl. 

Then the equations of the extremal lines Lt^^t2 can be rewritten in the form 

xi + 2tiAx3 + 2tiD = 
X2 + tlA'x3 + tjD' = 



(4.14) 



(4.15) 



(4.16) 



We need to work a bit to define functions A and D. We assume that our line 
intersects dU in a point C = C(*i)*2) on the "upper lid" X3 = M and in a point 
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^ = C(ti,t2) on the "side" X2 = x\. Then we have two pairs of equations ()4.16p 
asserting that our points are on the hne Lti,t2 

Ci + 2tiAM + 2tiD = (4.17) 

C2 + tjA' M + tlD' = (4.18) 

6 + 2tiA^3 + 2tiD = (4.19) 

+ tjA'^s + tjD' = (4.20) 



and two boundary conditions ()4.2p and ()4.3p 

^Cih + C?t2 = C?6 (4.21) 
ylt? = Cf . (4.22) 

Thus, we have six equations with six unknown functions: Cii C2, Cii ?3) A, and 
D. Equation (j4.18p determines the function (2- If we take Ci from (j4.17p and plug 
into K22t\ we get 

A = 4.{AM + D)2 . (4.23) 
So, introducing a new function of t2 

a = 2{AM + D) (4.24) 
we can express both functions A and D in terms of a: ()4.23p yields 

A = (4.25) 
and directly from the definition of a ()4.24p we get 

D = ^a- Ma^ . (4.26) 

Note that deducing (j4.23p we divided both parts of (j4.22p over ti. This is a 
correct operation because the Bellman function B clearly depends on all variables 
Xi, and therefore its partial derivatives ti cannot be identically zero. By the way, 
since is not identically zero, so is a. Moreover, we can assert that a is not a 
constant function. Indeed, assuming a to be a constant, we have C2 = from (I4.18|) . 
whence Ci = due to Ci ^ C2 = 0. But (j4.17p can be written as Ci + cih = 0, and 
since neither ti nor a is not zero, we come to a contradiction. 

In result we reduce our system to the system of three equations with three 
unknown functions 0(^2), Ci(^i)^2), and C3(ti,i2)- Indeed, equations (|4.19p . (I4.20p . 
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and (I4.2ip can be rewritten as follows 

ei+tia(l-2a(M-e3)) =0 (4.27) 
e? + t?a'(^-2a(M-e3)) =0 (4.28) 

^6ii-C?(e3-t2) = 0. (4.29) 

To have possibility divide ()4.29p over ^i, we need to check that d = cannot be 
an appropriate solution. Assuming that identically = we get 1 — 2a(M— ^3) = 
from (j4.27p . which turns (j4.28p into —^tfa' = 0. But this is impossible because, 
ti 7^ and a is not a constant function. So, we can rewrite ()4.29p as 

6 = ^2 + ^. (4.30) 

In fact ^3 does not depend on ti. Indeed, if we introduce 

rj=l- 2a{M - 6) , (4.31) 

equations (|4.27p and (|4.28p turn into 

^1 + tiar] = (4.32) 
e? + t?a'(77-^) =0, (4.33) 

whence 

a'(r?-^)+aV = 0. (4.34) 

Since a does not depend on ti, the function r] (as a solution of this equation) 
also depends only on t2- Soon we shall show that X3 does not depend on t2 either, 
it is a constant function! 

But now we remove ^3 temporary from the play. From (j4.3ip we have 

6 = M - , (4.35) 

la 

and (|O0]) together with (021) yield 



The resulting equation is 



or 



6 = t2 - . (4.36) 
lar] 



t2 = M '- , 

2ar/ 2a 

t2-M= + ^ . (4.37) 

Zarj 
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In result, we have two equations ()4.34p and (14.370 for two unknown functions a and 
rj of the variable t2- Of course, now we could solve quadratic equation (I4.34p with 
respect to rj and plug the solution in (j4.37p trying to solve the resulting differential 
equation with respect to a. This is possible, but it leads to many complications. 
Much simpler is to differentiate (|4.37p 

^ avi'{rf — 1) — a'irf — + rj) 

and then to replace the denominator by using (j4.34p 

a'{2ri - 1) = ar]'{rf - 1) - a'{rf' - rf + rf) , 

or 

{ri^ - l){r} - l)a' = a{ri^ - l)rj' . (4.38) 

First we consider the possible constant solutions 7/(^2) = ±1. Solution = 1 is 
not suitable. Indeed, if r/ = 1, then from (I4.37P we get 

1 



2(t2 - M) ' 



whence 



9 1 1 o t2 - 2M , , 

^ = « = J^^^^ ' ^ = 2 « - = 4(^^3^ ' (^-^9) 

A' = ^ D' = - '^-^^ (440) 

2(t2 -Mf 4(t2 - M)3 ■ ^ ^ ^ 

Therefore, in this case, the extremal line Lt^^t2 (cf- (|4.16p has the form 

tixz t2-2M _ 

+ 2(i2-M)2 + *^2(t2-M)2 - ° 

_ tfx3 _ 2 t2-3M ^ 
2(f2-M)3 M(t2-M)3 

But the only point of this line with = M belongs to our domain $7, because for 
other points we have 

2_ tf(2x3 + t2-3M) tl{x3 + t2-2M)^ _ tf(x3-M)2 
4(t2-M)3 4(t2-M)4 4(t2-M)4 

Solution ?7 = — 1 is also impossible. Therefore, we come to 

(rj - l)a' = ari' , (4.41) 
which implies a = C ■ {ij — 1). Then (|4.35p yields 

^3(^2) =M + ^=:c. (4.42) 
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Recall that ^3(^1, t2) = c is the level on which the line Lt^^t2 intersects the boundary 
X2 = Xi, < x-^ < 1, i.e., the third coordinate of the point of intersection. So, 
m < c < M. By our assumptions these lines intersect the boundary 2:3 = M, 
therefore, to foliate the whole domain Vt (at least its interior) we need to have 
= 171. So, from now on we assume ^3 = c = m. Of course, we could consider any 
c > m and try to find an additional solution of the Monge-Ampere equation in the 
subdomain m < x-^ < c but we shall not do this. First of all, we do not need to 
look for another solution. The better way of doing is to try now to check that the 
found function is just the Bellman function of the problem. Would we have some 
obstacles in proving we could continue our search. Another argument explaining 
why we do not need "to glue" our candidate from two parts is the fact that there 
is no concave solution in the subdomain m < X3 < c being continuous extension of 
the found solution. 

Recall that we have the following expression for B (see (I4.10p ) 

B = ^tixi + t2X2 . (4.43) 

Thus, to find an expression for B we need to find ti and t2 as functions of a point 
X running over our domain 
From (j4.36p we have 

t2 = m-\ , TT . (4.44) 

2a[l - 2a(M - m)] ^ ^ 

This gives us the desired a as a function of t2, but it is clear that a is more 
convenient parameter than t2- We shall express all other functions of t2 in terms 
of a and look for a as a function of x G 17. To this aim we return to the equations 
of the extremal lines (j4.16p rewriting them in terms of a. From (j4.44p we have 

dt2 1 - 4a(M - m) 



therefore 



and 



da 2a2[l -2a(M-m)]2 



2a^[l-2a(M-m)P 
~ 1 - 4a(M - m) ^ ' 

2 n / 4a3[l-2a(M-m)]2 

A = a%A' = 2aa' = } , (4.46) 

1 - 4a(M - m) ' ^ ' 



n 1 -n' ^ '(^ A m\ a^(l - 4aM)[l - 2a(M - m)]^ 

D = -a-Ma ,D =-a{l-AaM) = ^ _ ^a{M - m) ' ^^'^^^ 
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Now equation ()4.16p of the line Lt^^t2 becomes 
Xl = -tia [1 - 2a(M - x^)] , 

2 2[l-2a(M-m)]2[l-4a(M-X3)] 

X2 — 
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(4.48) 



^ 1 - 4a(M - m) 

As we have known, the side surface {x2 = x^} is intersected by -Lti,t2 the points 
of the "bottom hd" {3:3 = m} 



^i = -tia [1 - 2a{M - m)] . 



(4.49) 



If we take ^1 for parameterizing the extremal lines, then their equations will be 
more symmetric. Indeed, ()4.49p implies 

6 



1 + 2a(M - m) ' 



(4.50) 



and therefore ()4.48p turns into 



_ 1 - 2a(M - X3) 
~ 1 - 2a{M - m) ' 

1 - 4a(M - xs) ^2 

^2 = T 4l • 



(4.51) 



1 - 4a(M - m) 

These equations immediately supply us with an expression for a, namely, a = a{x) 
is a root of the following cubic equation 



s :- 



X2 



1 - 2a{M - xa) 



n 2 



1 - 4o(M 



m) 



1 - 4a(M - X3) 



(4.52) 



1 - 2a{M - m) 

To determine, which of three possible roots of the equation (14.521) gives us the 
desired value of a, we investigate the above function s = s{a) defined by ()4.52p as 
a function of the parameter a, all other parameters assuming to be fixed. 

First of all we note that the extremal line (j4.5ip intersects the plane {X3 = M} 
at the point 

6 



C = (Ci,C2,M), Ci = -_2,(^_^), 

Since C,2 ^ 0, we have the first restriction for a: 

1 



C2 



1 - 4a(M - m) 



a < 



(4.53) 



(4.54) 



4(M - m) 

The behavior of s{a) on the semi-axis (|4.54p is the same for all values of other 
parameters: on the negative half- line it monotonously increases from till 1, 



and on the interval [0 



' 4(M-m)- 



it monotonously decreases from 1 to 0. Therefore, 
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for a given s, < s < 1, we have one or two sohitions of ()4.52|) satisfying ()4.54p : 
for aU s there exists a positive solution a = and for s > "^"^ there exists a 
negative solution a = a~ as well. We shall see that our solution is a = a"*", but 
the solution a = a~ is not meaningless: unexpectedly we get a solution of another 
extremal problem, namely (jl.lOp . 

To complete step II of constructing a Bellman function candidate we need to 
write down formula (jl.Sp or (jl.lip . Everything is ready to do this: we have ex- 
pression (|4.43p for B, where we need to substitute (|4.44p for t2 and to take ti from 
the first hue of (|4.48p . So, we get 

B(x:m,M) = ^ + ( m -\ J\ X2 (4.55) 

^' ' ^ 2a[l-2a(M-X3)] V 2a[l - 2a{M - m)] J ^ ' 

(xs - m)x2 , . 

= K h mx2 , (4.56) 

[1 - 2a(M - m)] [l - 4a(M - xg)] 

where a = a"*" is the positive root (4(M — m)a+ < 1) of the cubic equation (j4.52p 
for Theorem [3l and a = a~ is the negative one for Theorem [H 

□ 



5. Foliation according to B. 

We already saw that lines ^ given by equations (j4.5ip and in the parametriza- 
tion of the family of the extremal lines it is more convenient to use parameter a 
instead of t2 and instead of ti: 

For a given we have a "fan" of extremal lines starting at an arbitrary point 
^ = iCij^ij''^) on the "edge" {x: X2 = xf, X3 = m}: 

-^iT = {f-* ■ - < a < 0} , 

where is the fan of extremal lines of -Bmax and F^^ is the fan of extremal lines 
of Bmm ■ 
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C2 



Let us find the trace of our fan on the boundary = M. We remember that 
the point of intersection ( is given by (I4.53|) : 

V = 

1 - 2a{M - m) ' 

1 - 4a(M - m) ' 
This is a parametric equations of the hyperbola C^^ 

(2 = ^^, (5.1) 

tangent to the parabola ^2 = Ci at the point (^i,^^)- 

The intersection of the fan with the plane X3 = M is the piece of this 
hyperbola between the points ^1 and 2^i, for the fan it is the piece of this 
hyperbola between ^1 and zero. For ,^1 = it is the axis = {(0, C2)}- 

Everything is on the "upper lid" {xs = M,xf < X2} of Q, and the curves C^^, 
as well as the curves C~^, foliate this "upper lid". Lines ia,^i G connect the 
points of to the point {^i,(,i,m) on the boundary of 0, and each of two sets of 
lines : — 00 < ^1 < +00} foliates some subdomain of Q. The lines of the fan 
Ft foliate the whole $7, whereas Fr foliate the subdomain 



{x = {xuX2,X3}:^-^<^<l, xs<M}, (5.2) 
^ M — m X2 

where the function -Bmin is defined. Indeed, we have a solution a = of (j4.52p 
for arbitrary point of 0,, while a solution a = a~ exists if and only if the point x is 
from ([O]) . 

6. Lower boundary of Q. 



Notice that we found our function B only in {xf < X2, m < x^ < M} for -Bmax 

2 

and in {xf < X2, M - {M - m)^ < X3 < M} for B 

iiiin- However both these 

functions are continuous in the closed domains. For Smin clearly the limit is 171x2 
as X3 ^ M — (M — fn)-^, what corresponds to a ^ —00. For -Bmax formula (j4.56p 



has indeterminancy g when 3:3 — > m, what corresponds a — > ^]^_^^ ■ But we can 
easily pass to the limit in ()4.55p and obtain by continuity the values of the function 
Bmaxixi, X2,m;m, M) on the "lower lid". Namely, 

Bma.xixi,X2,m;m, M) = 4(M - m)(x2 - x\) + mx2 ■ (6.1) 
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The norm of embedding becomes apparent. In fact, tiie best constant in the 
inequaUty 

Bmax 1^ C X2 

can be read from (j6.ip . and it is indeed 4(M — m). This can be easily proved, but 
notice that it is not even explicit in the formula for Bmax^- 

7. Reducing parameters M and m. 

Before we start to prove Theorems [3] and [H we would like to show that it is 
sufficient to prove them for m = and M = 1. On the one hand, we have 

B{xi,X2,X3; m, M) = B(xi, X2, — ; 0, 1) + mx2. (7.1) 

M — m 

This is direct consequence of definition. Indeed, while the homogeneity condi- 
tion (j4.4p follows from the definition if to compare a set of test functions {(/>: {(j)) j = 
{4'^) I = ^2} and the set {(f) = t(fi] with the same measure and the same set 
of point masses aj, relation (j7.ip follows from considering the same set of test 
functions (f) but comparing with a "renormalized" measure jl{J) = ^^AZ-m*^^ ^'^'^ 
point masses aj = jji^^ then we shall have m = 0, M = 1, and X3 = 

Notice that function B(x) from p.5p is defined differently than B(7£;7^(j;; 0, 1). 
Its definition does not allow for measure /i on the boundary. However, these func- 
tions turned out to be equal, extremal /i must be zero. It is not clear how to see 
this immediately from the definitions (jl.Sp and (jll.ip of B(x) and B(7£;7^(x; 0, 1) 
correspondingly. However, since the supremum in the definition of B(x) is taken 
over the smaller set of test measures [fj, = 0), we have the inequality B(x) < 
^cet{x]^, 1). So, the proof of Theorems [5] and [3] will consist in proving two in- 
equalities: Bcet{x] 0, 1) < B{x; 0, 1) and B(a;) > B{x] 0, 1). 

Let us introduce the notations: 

B,{x) := B{x; 1 - c, 1) , B,{x) := BceHx; 1 - c, 1) c G (0, 1] . (7.2) 

Obviously we can rescale everything and consider only Be and Be. Moreover, it 
is enough to consider only c = 1. 

In fact, our Bellman function candidate B{x;m, M) clearly satisfies relation (|7.ip . 
therefore, to prove BcETix',Tn, M) = B{x;m, M) it is sufficient to check that 

Bcet{x;0,1) =B{x;0,1) 
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for all suitable arguments x. We will first prove that 

B{x) > B(x;0, 1). 

Then we prove 

B <B. 

Then the obvious inequality 

B < B(x;0, 1) 

finishes the proof. 



8. Bcix) > Bc(x). Concavity. 

In what follows symbol B{x) stands always for B{x;0,l) from ()1.6p (which is 
the same as Bc{x) with c = 1). 

Theorem 5. Function Be is concave in the domain f^c- 

Proof. It is enough to prove this for c = 1. Then rescaling proves the rest. In what 
follows a is always the unique root of (j4.52p (with M = l,m = 0) lying in [0, |). 
In previous sections we calculated (c = 1): 

dB xi 

d^i^^^^ ~ a[l - 2a{l - xs)] ^ ' 

dx^~^^~ 2a(l - 2ca) ^^'^^ 
dB _ _ xj 

a^-'^- [l-2a(l-X3)]- ^^-^^ 
We can now compute Hessian d? B, where (we use the notation s = x\/x2) 

B{x) := M{x) = 

1 I l-4a(l-3;3) 2x1 da l-4a xi da 2x1 _ xy[l-Aa{l-X3)] _da_ 

' a[l~2a{l-xs)\ a2[l-2a(l-a;3)]2 ' X2 ' ds ^ a^(l-2a)^ ' X2 ' 9s ' [l-2a{l-x-i)]'^ a'^[l-2a{l-x:i)]'-^ ' dxi 

l-4a(l-X3) xf da l-4a a;f Qa l-4a da 

a^[l-2a{l-a;3)]^ ' 1^ ' ds ^ 2a^(l-2a)^ ' ^ ' ' 2a^{l-2a)^ ' 9x3 

2x1 _|_ 1-X3 8x£ aa 1-X3 4x4 ^ 4xfa 4xf(l-X3) 

[l-2a(l-X3)]2 "I" [l-2a(l-X3)]a ' X2 [l-2a(l-X3)]» ' ^ ' 5x3 [l-2a(l-X3)]^ ~'~ [l-2a(l-X3)]3 ' 8x3 

The element M12 actually is equal to M21, M13 is equal to M31 by ()4.52p . The 
reader may try to prove that this matrix is non-positive for every x € hy direct 
calculation. We prefer an oblique way of doing that. Put 
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N{x) :-- 



+ 



1 - 4a(l - xs) 



2x2 



a[l-2a(l-X3)] a2[l - 2a(l - X3)]2 X2 



da 
ds 



1 - 4a 



Xi 



l-4a(l-X3) 
'a2[l-2a(l-X3)]2 



9a 
9s 



a2(l - 2a)2 X2 
xf 



1 -4a 



2a2(l -2a)2 



L{x) :-- 



2x2 



l-4a(l-X3) 

a2[l-2o(l-X3)]2 X2, 

1 — 4a(l — X3) X' 
~a2[l-2a(l-X3)]2 ■ ^ 



9a 

9s ' 
9a 

9s 



1 - 4a 



a2(l - 2a)2 
1 -4a 

2a2(l-2a)2 



Xi 
^2 



da 

ds 
da 

ds 



Let us prove that 



Mii<0, L<0, N <0. 



We know that function s ^ a(s) decreases from j to when s goes from to 1. 
Therefore, a' > 0. Also 1 — 4a(l — X3) = (1 — 4a)+4ax3 > because a < ^, X3 > 0, 
and a > 0. Therefore, Ln < 0, A'^22 = -^^22 < 0, Nn < for every x G iln {0 < X3}. 
In particular, Mn < 0. 

On the other hand, it is immediate to see that L11L22 — -^^12^21 vanishes identi- 
cally. Then we can see right away that < as well, but more than that 



N <0, det iV(x) / Vx G 17 n {0 < X3} n {xi / 0} . 



■4) 



In fact, it follows because N = L + diag{- 



a[l-2a(l-X3)] 



< 



, 0} if we notice that 



a[l -2a(l -X3)] 
and L22 < 0. 

We would like to recall the reader that we obtained B of (|4.55p - (|4.56p just by 
solving the Monge-Ampere equation det (PB = in 17. So, of course, det M(x) = 0. 
However, we propose to the reader to check this as follows. 

Fix a point x G 11,0 < X3, and consider the line Lt^^t2 passing through x. 
Its directional vector was computed, it is (see (|4.13p . (|4.15p . (j4.46p . and (j4.48p 
combined): 

4a(l - 2a)2xf 



d{x) 



2axi 



1 



(8.5) 



-l-2a(l-X3) ' (1 -4a)[l - 2a(l -X3)]2 

Actually we huilt B from the condition that (fiB annihilate the vector field 
(i(x),x G H. Also it is easy to see that M{x)d{x) = from the direct calculation. 



da 

ds 
da 

ds 
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We saw that for every x G Q D {0 < x^} we have Mn < 0, M11M22 - < 
0, det M(x) = 0. By a well-known fact from linear algebra we conclude that B = 
M{x) is negatively defined for all such x. Concavity of B is fully proved. 

□ 

Theorem 6. In 0, we have §^ix) > xf. The same is true for Be in $7c for 
CG (0,1]. 

Proof. The second claim follows from the first by rescaling. The first claim is 

obvious from (18.10 : — — = --k > x? because < 1 — 2a(l — X3) < 1. 

dx3 [1 - 2a(l - xsjj^ 

We already noted the left inequality, it follows from a > and a < ^ < ^. The 
right inequality is just a > 0. 

□ 

Theorem 7. Let c £ (0, 1]. In Qc have Be < -Be- 

Proof. Again it is enough to consider only the case c = 1. Let us combine Theorem[5] 
and Theorem [6] to obtain the following main inequality: 



B{xi,X2,X3) - ^[B{xf,x^,x^)+B{x-^ ,X2,x^)]> xl[x3 - ^{x^ )] , (8.6) 
Vx, x~^,x~ G such that xi = ^{xf + x^) , X2 = ^{^2 + ^2 ) ' (^■'^) 

X3 = \{xt + x^) . (8.8) 
In fact, this is easy. Put Xi{t) = ^[xf{l + t) + x^ {1 — t)] , i = 1, 2, and 



X3(t) 



-txj +X3(l + t), if t£ [-1,0], 
X3(l-t)+t3;+, if t G [0,1]. 



Set 



h{t) := B{xiit),X2it), xsit)) , t G [-1, 1] . 
Then the main inequality above transforms into 



b{0) - ^(6(1) + H-l) > xi(0)2[x3(0) - 1(X3(1) + X3(-l)] . (8.9) 
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To prove (18. 9|) let us notice that X3(t) is concave, and negative measure x'^ is the 
fohowing x'l{t) = — 2[x3(0) — |(x3(l) + X3(— 1)] • 5q. Then, of course, h is concave 
(see Theorems [5] and [6]) and measure 

h"(t) = {<fBx'{t),x'{t))dt + ^{x{t))xl{t) = {(fB{x{t))x'{t),x'{t))dt 

ax3 

-2^(X(0))[X3(0) - ^(X3(l) + X3(-l)]5o . 

The following formula finishes the proof: 

h{^)-\{b{l)+h{-l) = -\ j' {I - \t\)h" {t) . 
In fact, combining the last two formulae we obtain (we use Theorems [5] and [6] again) 



6(0) - \{h{l) + h{-l) = \ |\l - m-d^B{x(t))x'{t),x\t))di 

+ ^(X(0))[X3(0) - i(x3(l) + X3(-l)] 
>Xi(0)2[x3(0)-i(x3(l)+X3(-l)], 

which is desired (|8.9p . 

Now the main inequality ()8.6p and the convexity of the domain will 
allow us to finish the proof of the theorem. We need a simple lemma about sub- 
harmonic functions on graphs. 

We consider the dyadic tree T, whose vertices are denoted by Po-, where a is 
the word formed by it, the empty word do is the root of the tree, and \a\ is the 
length of the word. Function on T is called superharmonic if its discrete Laplacian 
is non-negative 

E T A/(a) = /(a) - ^(/(a+) + /(a-)) > . 

We can associate the boundary of the tree with segment [0, 1] in the following 
sense. For Lebesgue almost every point x G [0, 1] we have a unique branch h(x) = 
(do, cTi(x), cr„(x), ...) of T associated with it: just consider the dyadic form of 
X and its n-th digit encodes whether to branch to + or — side on the n-th stage. 
Given / on the tree we put 

F{x) = liminf /(c7„(x)) . 

n— >oo 
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Here is a lemma which deserves to be cahed Green's formula for the dyadic 
tree 

Lemma 8. Given a positive finite superharmonic function f on T we get 



^0 "T^- 



crer 

Proof. Obvious. □ 

Now fix X = (xi, X2, X3) E O and fix any function (j) on the interval I = [0, 1] (we 
recall that nothing depends on the original interval) with 

(</.)/ = XI, {cj)^)j = X2, (8.10) 

a non-negative measure fi on /, such that < /i(J) < | J| for all dyadic subintervals 
J of /, and any collection of non-negative numbers {oj} jei:){7) such that 

JfF + Wl ^ "^ = '^=^^^°'^]' + ^ ^ aj e[0,l]W€D{I). (8.11) 

' ' ' ' J£D{I) ' ' ' ' JeD{£) 

We immediately see that /u = wdx, where <w < 1. 

Intervals of D(I) and dyadic tree T are in natural one to one correspondence. 
We call /o- the interval corresponding to vertex a, I corresponds to (Tq- Consider 
the following function. Take B from (|4.55|) (with M = l,m = 0) and put 



M(a):=^(/.(/.)+ Yl 



and 



Then of course 



f{a) ■.= B{{<t>)i^,{r)i.,M{a)). 



fiao) = Bixi,X2,X3). (8.12) 

By the main inequality (j8.6p it is a superharmonic function on the tree T: more 
than that, we can estimate its discrete Laplacian from below. By this same (|8.6p 

Af{a)>{c^)lAMia). 

It is immediately seen that 



AM{a) = ajjm>0. 
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Therefore, combining the last two inequalities we get 



l-^O-Q I 



< A/(a) . 



One can observe that B{xi, X2, x^) > X2X3. Then 
/(a) > (02)/.M(a) > ■ 



\Ia\ 

Then for almost every x G [0, 1] we denote by cr(x) the branch landing at it and 
by (Tnix) the n-th vertex on this branch. Then for a. e. x 

liminf f{an{x)) > lim {4>^)i^ ■ = <f)'^{x)w{x) . 

n— »oo n— >oo \Ia\ 

Recall that I = I^q = [0, 1]. It is time to use Lemma [H the last inequalities, and 
Fatou's lemma (which is J liminf„ /„ < liminf„ J /„ for a sequence of non-negative 
functions /„): 

4>'^{x)w{x)dx + -^y^icpfj ai^ < B{xi,X2,X3) . (8.13) 



We used here our lemma and (j8.12p . 

But inequality (I8l3]) proves that B(a;;0, 1) < B(x;0,l). In fact, in the left 
hand side of (j8.13p we have arbitrary satisfying ()8.10p and arbitrary numbers 
a. satisfying ()8.1ip . Function B(xi,X2,X3) by definition is the supremum of such 
sums over all such functions and collections of numbers. Therefore, inequality 
Bc{x) > Be (a;) is completely proved and so is the theorem. 

□ 

9. Bc{x) < Bc(x). Extremal sequences. 
Theorem 9. Let c G (0, 1]. In {x^ < X2 , 1 — c < X3 < 1 we have Be > B^. 

Proof. Below /q = [0, 1]. We consider first only c = 1 and x = {xi,X2, 1) on the 
"upper lid" of J7. We want to fix a large integer n and to construct function 
and sequence {a"}^g£)(/Q) in such a way that the sum 

{(pn)ja'^ > B{xi,X2,l) - ^n, £n ^ ■ 

e&D{io) 

where 

{Mio = xi (9.1) 
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and 

i^Dlo = ^2 . (9.2) 

Fix a large integer n and split Iq = [0, 1] into the union of n + 1 intervals 1^ = 
j2-fc^ 2-^+1], 1 < A;<n, and J= [0,2^"]. MapsmA:(x) = 2'=j;-l map 4 , 1 < A; < 
n onto Iq. Consider preimages of Ig, s = 1, n, under m^, k = 1, n. We obtain 
Iks- Then o maps I^s onto /q. We also consider J^,, A; = 1, ...,n, preimages 
of J under m^. We iterate this procedure obtaining /fcifc2...fcm' Jkik2...km,-i- Put 
a^^^ := 2-'= • 2-" = I Jfcl , < A: < n. We put ^^^^ = | Jfcifc2...fc„|. For all other 
dyadic intervals we put a/ = 0. We call J the generation, and all Jkik2...km the 
m-th generation. Obviously sum of lengths of all generations of J's is equal to 1. 
There for we have proved 

^ ae = l. 

££D(Io) 

Moreover, it is easy to see that all J's are disjoint (except may be the end-points) 
so 

ai< \I\ V/ G D{Io) . 

ieDiio),ici 
Now we are going to define 



Cyi 1 *^ ^ J 

(/)„(mfc(x)) X € Ik k = 2, ...,n (9.3) 
d„(/)„(mi(x)) X e h. 



Here c„ and dn are constants which we will define now. First of all notice that this 
recursive definition of (pn really defines it a.e. (as long as Cn,dn are prescribed) 
just because we have it already defined on the generation of J's (that is on J 
itself). But then it is defined on m^-pre-images of J (so on the first generation of 
J's), but then it is defined on pre-images of pre-images (second generation of J's), 
et cetera... But The union of all generations of J's gives us Iq up to a set of zero 
Lebesgue measure. 

We define c„ , dn from requirements (|9.ip , (j9.2p : 



1 , 1 1, 1 

— CnXi + (1 - — - -)Xi + -dnXi = Xi 
— 0^X2 + (1 - — - -)X2 + -d„X2 = X2 . 
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Then we get (using our notations s = x\/x2 and choosing proper root of qua- 
dratic equation) 



Lemma 10. We have 

Hx)CnXJk,...k,„ < Cn(</'n)/,^...fe„ , Vx G Jfci...fc„ , (9.5) 

where empty sequence ki...km corresponds to Iq and J correspondingly. 
Proof. Let us start with the empty sequence. Then we should check that 

H^) < Cn{(t>n)lo , Vx G J. 

But the average over Iq is xi, and (j){x) = CnXi on J by definition. The rest 
of the lemma follows from the self-similarity of in fact, it is either the same 
function (pn "pre-shrunk" to a smaller interval, or it is a fixed multiple of this. So 
we are done. □ 

Let us now square (|9.5|) and integrate it. Then sum over all J's of all generations. 
We will get from ([O]) 

2„n ^ ^2 2 



esDiio) ^{l/s-^l/s^-l/sr 



But the first term is 

2 

XT 



(71^-71/^)2 



.?(v/T7I + yT77^)^ 



We finally get 

{Kfa'l > {V^+ \jx2- xlf - En . (9.6) 

teD(7o) 

Let us now look at formula ()4.56p with X3 = 1 and plug a from (j4.52p with 

2 

X3 = 1. Readily, (j4.52p with X3 = 1 becomes (^l^2a)'^ ~ ^2' -^'^o™ here and from 

negativity of a we obtain 



X2 - x\ 



= 1 + w 

1 + 2a \ X2 
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But formula ()4.56p with X3 = 1 (and c = 1) gives 




Using (j9.6p we finally we get 

^ {<pnfa'l>B{xi,X2,l)-en, (9.7) 

i€D{Io) 

and our theorem is proved for X3 = 1, c = 1. 

If c = 1 but < X3 < 1, then point x = {xi,X2,X3) G Q lies on a certain 
line I/ti.t2- Function B is linear on this line. So combining constant function and 
function (pn we will obtain the general inequality. 

For c G (0, 1) the proof is actually exactly the same. □ 

9.1. Geometry of foliation: explanation of the choice of 0„ in ()9.3p . This 
section contains a heuristic explanation why function 4>n was chosen in this form 
in (I93|) . 

The main point here is a certain geometric observation concerning the foliation 
of the domain 0, by lines ititj • Consider the foliation of the upper lid {X^ = 1} n O 
by parabolas Pa := {X2 = A Xf, = 1} , A > 1. Then we have the following 
easy geometric observation. 

Theorem 11. Let Lt^^t2 intersects the upper lid at the point {ui,U2,l) and let 
(ni,U2, 1) G Pa- Then the projection of the line Lt^t2 onto the upper lid is tangent 
to Pa at point (tii,U2,l). 

Proof. Fix ti,t2, which fixes a = 0(^2) and consider formulae (|4.48p that give us 
Lt^t2- Then we can compute the slope of the projection of Ltj^t2 on plane (Xi, X2). 
In fact from (|4.48p 

dXi _ 2 9X2 _ Atja^jl - 2af 

Therefore, the slope of this projection is the ratio of these quantities: 

2tia(l - 2af 

l-4a • ^^-^^ 

On the other hand the point of intersection of Lt^t2 wit the plane {X3 = 1} is 
(again see (I4.48j) ) 

U = (ui, U2, 1) , where ui = —tia , U2 = tf a^(l — 2a)^ . 



36 VASILY VASYUNIN AND ALEXANDER VOLBERG 

Therefore, point U lies on parabola ^3 = 1,^2 = ^^i-ta ' ^i' hence, the slope 
of this parabola at U is 

2(1 - 2af , , 

X__L.u.. (9.9) 

Now plug u\ = —t\a into (|9.9p to see that it is equal to the formula for slope in 
(|9.8p . We are done. 

□ 



Let us see what influence this theorem has on the choice of (^n- 

Given L = Lt^t2 intersecting the upper lid at u = (ni,M2,l). To construct 4> 
corresponding to u we make an infinitesimally small jumps u + ee, u — ee along 
L (here e is a unit vector parallel to L and with positive third coordinate). If we 
would know how to build the extremal functions (f)± for u it ee we would just put 
the one for u + ee on the right half /+ of / and the one for u — ee on the left half I_ 
of /. The one for u + ee is definitely unknown. But the one for u — ee is "known" 
in the sense that we can jump from it to n — 2ee = u — ee — ee and u = -u — ee + ee. 
Then we can restart the procedure again. Thus, the explanation for the self-similar 
structure of and for the second line of (|9.3|) . 

Notice that after several jumps along L we will find ourselves at the end-point 
u — nee of L where = 0, X2 = Xf (especially if we choose dyadic number e). At 
this moment WG know thctt cxtrGiricil function is just constant JCi — : Cn Ml. This is 
the explanation for the first line in the definition of 

To explain dn let us recall that we do not yet know how to build the extremal 
function for the point u + ee. Moreover, we cannot do that. This point is outside 
of fl. So here is heuristics. Instead of jumping into u + ee let us jump to point 

= projection of u + ee onto the upper lid. We know by Theorem [11] that 
lies on the tangent line to parabola Pa- X2 = Ax\ at point u. So, in a sense, we 
can think that Pe lies on P4 (because e is very small). Or rather we can choose 
U =pe + 0(e2) such that U G Pa- Let U = {Ui,U2,l)- 

Now we need to understand the following: suppose we know an almost extremal 
function for u (our future (pn)- Do we know it for U then? The answer is yes, be- 
cause if we know the extremal function for one point of parabola than the extremal 
function for another point of the same parabola is acquired just by multiplication 
of the first function on a suitable constant. 

This is the explanation for the third line of (|9.3p . 
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10. Step I for Carleson embedding theorems in with p > 1 

Now, we briefly consider how to find the Bellman function for the Carleson 
embedding operator acting on arbitrary L^. We shell follow the same scheme and 
we shall start with an evident definition of the Bellman function. 



{<p)j=Xi, {\4>\P)j=X2, jJi (/^(^) + X] "^)=^3, 



JeD(i) 

m\J\ < fi{J) + ae<M\J\ VJ G L>(I)} , (10.1) 



Bmin(xi,X2,X3;m,M) :=inf I / \(j){t)\P dfiit) + ^ \ {(P) jl^aj : 

' J€D(I) 

{4>)i=Xl, {\(I)\P) J = X2, jJi (/^(^) + X] "^)=^3, 

' ' J€D{I) 

m\J\ < fi{J) + ai<M\J\ VJgZ)(I)}. (10.2) 

i€D{J) 

Functions B are defined in 

:= {x = {xi,X2,X3) : x^ < X2, m < xs < M} . (10.3) 

As before, the Bellman functions B do not depend on the interval I. The differ- 
ence is that now we consider non-negative test functions only. 

From the definition of function B above we can immediately see that it satisfies 
the inequality (just concavity or convexity) 

Bmax(2;) - ^(Bmax(a;") +Bmax(2;+)) > 0, Vx± G , X = + x') . 



Bmm(2;) - ^(Bmm(x ) +Bmin(2;+)) < 0, ^X^ £ Q , X=^(x+ + x). 



Note that Q (see (|10.3p ) is convex and for all x^ £ ft we have x G 0. 
The boundary conditions also follow from the definition: 



(10.4) 

)• 

(10.5) 
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B(xi, |xi|^,X3) = \x1\Px3 . 



— {xi,X2,M;m,M) = \xi\p . 
0x3 



(10.6) 
(10.7) 



It is equally easy to see that homogeneity condition holds 



B{tXi,tPx2,X3) =tPB{xi,X2,X3), t > 0, . (10.8) 

As before, we are looking for a function B satisfying the following degeneration 
condition: 



det{(fB) = det 



Bx2X\ Bx2X2 Bx2Xz 



0. 



(10.9) 



\BxzXl Bx^X2 Bx 



3X3/ 



Lemma 12. There is a simple algorithm to find the function B{x;m, M) that 
solves Monge-Ampere equation (|10.9p in the domain (jlO.Sp with boundary condi- 
tion (|10.6p - (|10.7p and homogeneity condition (jlO.Sp . 



Proof. Again we look for a solution B of the form 

B{xi,X2,X3) =ti-Xi+t2-X2 + t3-X3+to 



(10.10) 



where t,; = -i^, i = 1,2,3. The integral curves of the vector field 6 such that 



dxi 



d'^B{x)Q{x) = are again segments of straight lines 

xidti + X2dt2 + 3:3^^3 + dtQ = . 



(10.11) 



Each line is given by fixing two free parameters, for example ti and t2- Then we can 
rewrite equation (jlO.lip defining lines Lj^^jj foliating our domain in the following 
form: 



Xl + 
X2 + 



dto 



dt3 

at2 



X3 + 



dto 
dt2 



Differentiating (jlO.Sp in t and setting t = 1 we get 

Xiti +pX2t2 = pB . 

Whence, 

{p - l)tiXi + pt3X3 + pto = 0. 



(10.12) 

(10.13) 
(10.14) 



MONGE-AMPERE EQUATION FOR CARLESON EMBEDDING THEOREMS 39 

If we subtract this equation from the first equation of (I10.12p multiphed by 
{p — l)ti, then we get 

{pts -{p- l)h ||)X3 + {pto -{p- l)tl = 

for every X3. If we introduce the dual exponent q = then the latter equation 
can be rewritten as follows 

/ dt^ ^ I dtQ \ 

[1^3 - h ^)X3 + [qto -ti—)=0. 



Hence 



qt3-h--^ = 
Oti 

dto 

qto-h^ = 0. 

Oti 



(10.15) 



(10.16) 



From these "PDE" we easily write down 

' t3 = A{t2)\h\'' 
to = D{t2)\tl\'' . 

Then the equations of the extremal lines L^^ can be rewritten in the form 

xi + q\ti\i-HiAx3 + q\ti\'i-HiD = 

(10.17) 

X2 + \tl\1A'x3 + \tl\'^D' = 

As before, we assume that our extremal line -Lti,t2 intersect the boundary in 
a point C = C(*i;*2) on the "upper lid" X3 = M and in a point ^ = ^(ti,t2) on the 
"side" X2 = Then we have two pairs of equations (|1U.17|) asserting that our 

points are on the line Lt-^^ 

Ci + q\ti\''~'^hAM + =0 (10.18) 

C2 + |ii|M'M+ |fi|'?L>' =0 (10.19) 

a + q\tirhiAC3 + q\ti\''-\D =0 (10.20) 

\(if + \ti\'^A'^3 + \ti\'D'=0 (10.21) 

and two boundary conditions (110.60 and ()10.7p 

-6*1 + 161^*2 =ieire3 (10.22) 
p 

A\ti\i =\Ci\^ . (10.23) 
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If we take Ci from (jlO.lSp and plug into (110.230 we get 

A = qP\AM + D\P. (10.24) 

So, introducing a new function of t2 

a = q{AM + D) (10.25) 

we can express both functions A and D in terms of a: (jl0.24p yields 

A = \a\P (10.26) 

and directly from the definition of a ()10.25p we get 

D = -a-M\a\P. (10.27) 
Q 

In result equations (|10.20p . (110. 2ip . and (I10.22P can be rewritten as follows 

6 + |ti|^"^tia(l - qa\a\P-^{M - ^3)) =0 (10.28) 
leir + - qa\ar\M - 6)) =0 (10.29) 

-6*1 - I6r(e3-t2) =0. (10.30) 

p 

As before, we introduce the following new function 

rj = I - qa\ar\M - ^s) , (10.31) 
and then equations ()10.28p and (|10.29p turn into 

+ \ti\'^-hiar, =0 (10.32) 

|6r + ^|ti|V(r?-^) =0, (10.33) 

whence 

a'(??--) + -|ari??r = 0. (10.34) 
q p 

Since a does not depend on ti, the function rj (as a solution of this equation) 
also depends only on t2- 

On the one side, from (jl0.3ip we have 

6 = M - ^ , (10.35) 

qa\a\P ^ 

on the other side (110. 30p together with (I10.32P yield 

e3 = t2 + 4^ = t2-^^. (10.36) 
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So, excluding X3 we get 

ari _^ l-V 



p\a\P\rj\P qa\a\P~'^ 

Differentiating in t2 the latter equation we obtain 

^ a'rj + rj'a rj' V ~ ^ (1 ~ 

g|a|P|r/|P qa\a\P~'^ q \a\P 

and then we multiply this equality by the common denominator and replace it 
using (jl0.34p by —pa'{ri — |): 

—pa'iri ) + a'rj + Vj'a = ??'a|r/|^ + (p — 1)(1 — r/)a'|r/|^ , 

or 

(p - l)(|r/|P -!)(??- l)a = a{\i]\P - 1)??' • (10.37) 
Thus, we come to the following equation 

-!)(??- l)a' = a??', (10.38) 

which implies r] = \ + C\a\P . Then pu.35|) yields 

^3(^2) = M + ^^fl = M + - sign a =: c. (10.39) 

The latter expression is a constant, because sign a can take the only value due to 
(^3 < M. As before, the natural choice of the constant c is c = m. 
Recall that we have the following expression for B (see (jl0.13p ) 

B = -tixi + t2X2 . (10.40) 
P 

Thus, to find an expression for B we need to find ti and t2 as functions of a point 
X running over our domain Q,. 
From (|10.3ip we have 

77 = 1 - qa\a\^''^{M - m) , 

hence ()10.36p yields 

ari a — qiM — m)\a\P , , , 

t2=m + = m+ ^ \r, - (10.41 

p\ar]\P p\a - q{M - m)\a\P\^ 

This gives us the desired a as a function of t2, but it is clear that o is more 
convenient parameter than t2- We shall express all other functions of t2 in terms 
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of a and look for a as a function of x G 0. To this aim we return to the equations 
of the extremal lines (I10.17P rewriting them in terms of a. From ()10.41|) we have 

1 d , ^ 1 - pq{M 



da 



q\ar]\P da 



q\a - q[ 



therefore 



and 



{M-m)\a\p\^ 

qa\a — q{M — m)\a\P\^ 
a — pq{M — m)\a\P 



A 



D = -a 

q 



■ p\a\^ ^aa' 



-pq\a\' 



\a - q{M - m)\a\P\^ 
a — pq{M — m)\a\P 



M\a\P ,D' = -{l-pqM\c 



IP-2, 



)a' 



(a - pqM\a\P) \a - q{M - m)\a\P\ 
a — pq{M — ■m)\a\P 
Now equation (jlO.lTp of the line Lt-^^,t2 becomes 



t- 



[a - q{M - xs)\a\P) , 



X2 = \tl\' 



\a 



q{M - m)\a\P\^{a - pq{M - x^)\a\P) 



(10.42) 



a — pq{M — m)\a\P 

To get a more symmetrical expression, as before, we take a more geometrical pa- 
rameter ^1 instead of ii, where is the first coordinate of ^ = (^i, I'^il^j'n), the 
point where the side surface {x2 = is intersected by our extremal line. In 

term of a and equations ()10.42p turns into 



Xi 



X2 



a — q(M — X3)|a|' 
a — q{M — m)|a|i 

a — pq{M — 2;3)|a|- 



■6, 



16 r. 



(10.43) 



a — pq{M — m)\a\P 

These equations immediately supply us with an expression for a, namely and 
a = a{x) = is a root of the following equation 



s :- 



\xi\ 

X2 



q{M-x^)\a\ 



a — pq{M — 'm)\a\P 
a — pq{M — X3)|a|P 



(10.44) 



a - q{M - m)\a\P 

To determine which of possible roots of the equation (jl0.44p gives us the desired 
value of a, we investigate the above function s = s{a) defined by ()10.44p as a 
function of the parameter a, all other parameters assuming to be fixed. 
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First of all we note that the extremal line (I10.43P intersects the plane {xs = M} 
at the point 

c = (Ci,C2,M), ci = 7^ — C2 = Tw^^n^- (10-45) 

a — pq[M — xs)\a\P a — pq{M — X3)\a\P 

Since (2 ^ 0, we have the first restriction for a: 

a < . (10.46) 



^pq{M — m 

The behavior of s{a) on the semi-axis ()10.46p is the same for all values of other 
parameters: on the negative half-line it monotonously increases from (^E^)^ ^ 
till 1, and on the interval [0, {_ pq(^M-m) )'^ ^] monotonously decreases from 1 to 
0. Therefore, for a given s, < s < 1, we have one or two solutions of ()10.44p 
satisfying ()10.46p : for all s there exists a positive solution a = and for s > 
( M-m ) ^ ^ there exists a negative solution a = a~ as well. We already know that 
the solution a = corresponds to the function Bmaxj and the solution a = a~ 
corresponds to the function Bmin. 

To write down the expression for B we use (110. 13p . then we express xi using 
the first equation of (I10.42p . and by means of the second equation of (jlO.421) we 
replace the expression of \ti\''. For t2 we use p0.4ip . In result we obtain 

B = -xih + X2t2 = --|ti|^(a - q{M - X3)\a\'P) + 2:2*2 = 
p p ' 

ia - q{M - x^)\a\P)(a - pq{M - m)\a\A / a-q{M-m)\a\P 

— — Xo~\~Xo\ Tfi-\ — — ■ 

p\a - q{M - m)\a\P^\a - pq(M - x-i)\a\P) V p\a - q{M - m)\aY^ 

a|aP(x3 — rn)x2 

-|- mx2 = 



\a - q{M - m)\a\P\^ {a - pq{M - X3)|a|P) 



{xs - m)x2 



, —TYT, T + rnx2 ■ 

\l - q{M - m)a|a|P-2|^(l - pq{M - X3)a\a\P-^) 

As it was already mentioned, here we choose the solution a = a"'' of (jl0.44p for 
^max and a = a~ for B^nin- 

Theorem 13. Consider 

B{xi,X2,X3;m,M) = - (xg - m)a;2 + mx2 . (10.47) 

[1 - 2a(M - m)] [l - 4a(M - X3)] 
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We shall denote this function by B^^x if cl = a^{x) is the unique positive solution 
of the equation 



\xiY 

X2 



a - q{M - X'i)\a\P 



^ a — pq{M — m)\a\P 
a-pq{M - xz)\a\P 



a - q{M - m)\a\P 

from the interval [0, ( pg(jj_m) ^1' ^^'^ ^^^'^ function will he denoted by -Bmin if 
a = a~{x) is the unique negative solution of the same equation. The domain of 

Bmax is 

Dom(Smax) = {x: X2> \xi\P, m < xs < M} , 
the domain of Smin is 

M — m \p-i. 



Dom(i3mm) = {x-- \xif < X2 < 



\xi\P, m<X3<M} 



Then 



Bmaxix;m,M) = < 



and 



Bmin{x;m,M) = < 



M - 

mx2 X3 = m . 

Srnin(x;m,M) M- (M-m)^ < X3 < M, 



mx2 



m < X3 < M - (M - m] 



□ 



The proof of this theorem follows almost verbatim the lines of the proof of our 
Theorems EJ H] above. So, we are not going to prove it here. Instead, we wish to 
show that our results encompass a slightly more general situation than it could 
have been thought. 



11. Trees. 

As in [M] we let {X, fio) be a nonatomic probability space. Two measurable 
subsets A, B oi X will be called almost disjoint if fio{Ari B) = 0. Then we give 
the following 

Definition. A set T is called a tree if the following conditions are satisfied 
• X £ T and for every / G T we have iJ-o{I) > 0. 
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• For every I £ T there corresponds a finite or countable subset C{I) C T 
containing at least two elements such that: 

a) the elements of C(/) are pairwise almost disjoint, 



It is clear that E has /io measure zero. 

An easy induction shows that each 7^ consists of pairwise almost disjoint sets 
whose union is X. Moreover, x £ X \ E, then for each m there exists exactly 
one Im in %n containing x. For every m > there exists a J £ %n-i such that 
Im £ C{J). Then of course J = Im-i- Hence the set A{x) = {I £T : x £ 1} forms 
a chain of nested elements of Tm,m = 0, 1, .... Inclusions are all strict. From this 
remark it clear that if /, J G T and I D J Ci {X \ E) is nonempty then / C J or 
J G I. In particular, for any I, J £ T we have either n J) = or one of them 
is contained in the other. We denote by T( J) all elements of T that are subsets of 
J. Also for any I £ T put 



Given any tree T consider the collection a = {a/}/g7- of non-negative numbers. 
Let be another measure on X. We can introduce the tree Bellman function: 



T = Um>o'7'm, where % = {X}, TT^+i = ^ie.r^C{I). 
We have limm^oo sup^g^^ /Uo(/) = 0. 



For any tree T we define its exceptional set E = E{T) as 



E = U/gT Uji,j2ec(/),Ji^J2 -^1 n J2 . 






miioiJ) < + ^ < M/xo(J) \/J £T 

t&r(j) 



}. (11.1) 



Independently of the tree the same Theorem [3] holds: 



Theorem 14. Consider 



B{xi,X2,X3;m,M) 



{xs - m)x2 



+ mx2 



(11.2) 



[1 - 2a(M - m)] ^ [l - 4a(M - X3)] 
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where a = a{x) is the unique solution of the cubic equation 

x\ _ [1 - 2a(M - xa)]^ 1 - 4a(M - m) 
X2 ~ |_ 1 - 2a(M - m) J 1 - 4a(M - xs) 



(11.3) 



on the interval [0, 



1 



Then 



A{M-m) 



< 



B{x;m,M) 



mx2 



X 



1 < X2, 



m < X3 < M 



x^ = m. 



We can introduce the tree Bellman function also for p ^ 2,1 < p < oo. Again 
the result for the trees will not depend on the tree, it coincides with dyadic result. 
The proofs of Bj^2^(x;m, M) < B{x;m, M) is exactly the same as before, the 
proof of ^Q^rp^x; m, M) > B{x;m, M) is not too difficult either if one uses the 
following lemma from [M] . 

Lemma 15. For every I € T and every a such that < a < 1 there exists a 
subfamily F{I) C T consisiting of pairwise almost disjoint subsets of I such that 
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